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ABSTRACT 

Big bang nucleosynthesis constraints on baryon isocurvature perturbations are de- 
termined. A simple model ignoring the effects of the scale of the perturbations is first 
reviewed. This model is then extended to test the claim that large amplitude per- 
turbations will collapse, forming compact objects and preventing their baryons from 
contributing to the observed baryon density. It is found that baryon isocurvature per- 
turbations are constrained to provide only a slight increase in the density of baryons 
in the universe over the standard homogeneous model. In particular it is found that 
models which rely on power laws and the random phase approximation for the power 
spectrum are incompatible with big bang nucleosynthesis unless an ad hoc, small scale 
cutoff is included. 



1 Introduction 



Big bang nucleosynthesis (BBN) has produced well studied predictions of the light element abun- 
dances (Yang et al. 1984; Walker et al. 1991, hereafter WSSOK; Smith, Kawano, & Malaney 1993; 
Kernan & Krauss 1994; Copi, Schramm, & Turner 1994). These predictions restrict the total bary- 
onic contribution to the critical density of the universe, to be < 0.1. There have been many 
attempts to get around this bound and to extend it to the theoretically preferred value tto = 1. A 
notable class of such attempts is inhomogeneous BBN (for a review see Malaney & Mathews 1993). 
These studies of small scale inhomogeneities in the neutron to proton ratio including hydrodynamic 
effects, diffusion, extended networks, and multizone calculations turn out, however, to provide no 
appreciable increase on the bound to Qb set by standard BBN (Kurki-Suonio et al., 1990; Mathews 
et al, 1990; Terasawa k Sato, 1991; Thomas et al. 1994; Jedamzik, Fuller, & Mathews, 1994). 

On another front, structure formation theories are being constrained by a rapidly growing body 
of observational data. From this the primeval isocurvature baryon (PIB) model has fared relatively 
well. The PIB model relies solely on baryons to make up the matter in the universe and isocurvature 
perturbations to generate the structure (Peebles 1987a,b; Cen, Ostriker, & Peebles 1993). One short 
coming of this theory is that it requires f^o = SIb = 0.1-0.2, above the upper bound on fi^ from 
BBN. Independently, it has been suggested that non-linear isocurvature fluctuations may allow a 
larger contribution by baryons than allowed for in the standard, homogeneous case (Hogan 1978; 
Hogan 1993). 

In this work we have looked at the effect of large scale isocurvature perturbations on nucleosyn- 
thesis. Our treatment follows that of Epstein & Petrosian (1975, hereafter EP) and Yang et al. 
(1984). In those efforts, it was assumed that the volume distribution for the nucleon abundance 
could be described by a gamma distribution in the baryon to photon ratio, r/. The abundances of 
the light elements can then be used to constrain the parameters of the nucleon abundance distribu- 
tion. Here we will update that analysis utilizing the most recent constraints available from the light 
elements including 7 Li. We will also consider additional forms for the nucleon abundance distribu- 
tion to include the log normal (Barrows & Morgan 1983) and the gaussian (Sale & Mathews 1986) 
distributions. Furthermore we have extended the analysis to include a distribution of power on 
different scales to allow for dense regions to form compact objects and hence not contribute their 
light elements to the observed abundances. 

For our model we assume the perturbations can be described by a power spectrum with random 
phases. We have no knowledge of the spatial distribution, rj(x), and instead specify the choice of 
the density probability distribution, firj). Recently Gnedin, Ostriker, & Rees (1994) have also 
considered the effects of baryon perturbations on nucleosynthesis. They chose a log-normal density 
distribution with two scales. We choose three different density distributions (including the log- 
normal distribution) defined on a single scale. We have not specifically chosen parameters to match 
those of the PIB model. Our results, where comparable, agree with theirs. Gnedin, Ostriker &: 
Rees (1994) also approached the problem from the opposite direction; they assumed a form for 
rf(x) and derive a density distribution f(rj). They have found some models that can circumvent 
our bounds at the expense of assuming correlated phases. 

An outline of the paper is as follows: in §2 we discuss the observational bounds on the light 
element abundances used in this paper, in particular, how they differ from the values found in 
WSSOK. In §3 we describe our model for the inhomogeneities. In §4 we present the results of our 
calculations. 
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2 Observational Limits 



Observational measurements of the light element abundances play the crucial role of constraining 
the standard big bang nucleosynthesis model as well as models of nucleosynthesis which include 
inhomogeneities. The process of extracting abundances from the measurements, in particular pri- 
mordial abundances, is a difficult task. An analysis of this process in the context of limits on BBN 
is given in WSSOK. The 95% confidence limit (2a) primordial abundances quoted in WSSOK are 

Y P = 0.23 ±0.01, 

D/H > 1.8 x 10~ 5 , 

(D + 3 He)/H < 1.0 x 10 -4 , 

7 Li/H = (1.2 ±0.2) x 10~ 10 . (1) 

Here Yp is the 4 He mass fraction. These limits restrict the present value of r\ to 2.8 < 7710 < 4.0 
(7710 = 10 10 r/). In what follows we will use the WSSOK values with slight modifications to Yp and 
Li as discussed below. 

It is noted in WSSOK that the upper limit on Yp, Yp < 0.24, may be uncertain by 0.005. More 
recently a number of high precision measurements of 4 He in extragalactic H 11 regions have been 
made (Pagel et al. 1992; Skillman et al. 1994a,b; Izotov et al. 1994). Olive and Steigman (1994) 
have performed a detailed statistical analysis of these new measurements and found the primordial 
helium value 

Yp = 0.232 ± 0.003 ± 0.005, (2) 

where the statistical error is listed first and the systematic error second. The 95% confidence range 
(including systematic errors) is 

0.221 < Y P < 0.243. (3) 

We will employ this range in our analysis. 

The 95% confidence limit quoted by WSSOK for 7 Li/H consists solely of the statistical errors in 
the measurements. Recently Thorburn (1994) has made detailed measurements on a large number 
of metal poor dwarf stars. Her analysis employed a different model of stellar atmospheres than the 
one used to derive the data compiled in WSSOK. This model produces higher effective temperatures 
and hence higher lithium abundances. Her data yield a higher mean 7 Li abundance 

7 Li/H = (1.8 ±0.1) x 10~ 10 (4) 

where the quoted error is again only the statistical uncertainty in the mean. The difference between 
the 7 Li abundance given in (|l|) and (H) is a good estimate for the size of the systematic errors 
involved in making a determination of the the primordial 7 Li abundance. For this work we will 
consider both this new upper limit and the WSSOK upper limit. 

In summary, we are using the primordial abundances of D and 3 He as found in WSSOK (|T|) 
and modifications of the WSSOK values of Yp (||) and 7 Li @ due to recent measurements with 
explicit consideration of systematic errors. The primordial abundance limits used throughout the 
rest of this work are 
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The lithium bound is particularly important for constraining density fluctuations. In standard 
homogeneous BBN lithium must be nears its minimum value of 7 Li/H ~ 10~ 10 in order to be 
concordant with D and D + 3 He. Obviously any density variation selects 7 Li values above the 
minimum, hence 7 Li tightly constrains the range of perturbations. 



3 Model of Density Fluctuations 

We begin with a simple model of inhomogeneities (EP). We assume that some unknown process 
generates a baryon to photon ratio rj(x) at each point in space in such a way that the fraction of 
regions with a given value r\ is governed by the distribution f{rj). We acknowledge our ignorance of 
the process that generates i](x) by assuming (instead of deriving) the form of f{rj). Given f(rj) each 
region has a constant baryon to photon ratio r/ throughout nucleosynthesis. The regions undergo 
standard BBN, then mix producing the observed abundances of light elements. The distribution 
of these regions is described in our model by the function f(r]). Given the distribution f(rj) the 
average mass fraction of a light element is 

Xi= J™ d-nrifWXiW J fj, (6) 

where Xi{rj) is the mass fraction of element i according to standard Big Bang nucleosynthesis in a 
region with a baryon to photon ratio of r\. The average value of r\ for the universe is given by 

fj= dq 77/(77), (7) 
J 

for f(rj) normalized, 



r 

Jo 



d V r,f(r,) = l. (8) 

/o 

Shortcomings of this simple model include the assumption of equal power in perturbations on 
all scales and the allowance of extremely dense regions to contribute to the observed light element 
abundances today. It has been pointed out (Rees 1984) that high amplitude perturbations with 
a mass larger than the Jean's mass at the time of recombination will form gravitationally bound 
objects. These objects prevent the baryons in them from mixing with other baryons in the universe. 
Hence these overdense regions would not contribute to the observed abundances. 

To determine a more realistic model, we consider isocurvature perturbations to the baryon to 
photon ratio. These perturbations are characterized by their power spectrum, (|5fc| 2 ). Given the 
power spectrum the average abundance is found by 



Xi = 



f<^(\6k\ 2 )fo Vc{k) dvvf(v)Xi(v) 



S^{\Sk\ 2 )S3 eW dnvm 



JfA^k)J^ k) d vv f( v ) 



Here 



&\k) = ^(\6 k \ 2 }, (10) 

rj c (k) is the cutoff in rj based on the Jeans mass at recombination, and /c max is imposed to insure 
that the integrals converge. Since the Jeans mass at recombination, Mj oc r/ -1 / 2 (Hogan 1978; 
Hogan 1993) and the mass inside a scale k, M k oc i]/k 3 , the cutoff in i] is 

rj c {k) oc k 2 = (5k 2 , (11) 



(9) 
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where j3 ~ 6 x 10 for k in Mpc . Using this relation and interchanging the order of integration 
in Xi (|) gives 

where r/ max = ^c(^max)- In Eq. ( |l2|) our two assumptions are manifest. We have imposed an upper 
limit on k max implying that there is no power in the perturbation spectrum on scales smaller than 
A m in= 27r/fe max and we have assumed that on sufficiently large scales, corresponding to > Mj 
or k < yry//3 gravitational collapse will prevent these regions from mixing the hence these regions 
to not contribute in an average element abundance. One should note however that the average 
value of 77, 77 is not constrained by 7/ max . The density distribution indeed includes regions with 77 > 
f?max, though they do not contribute to the quantities Xj. We will rewrite A; max in terms of r/ max 
in the rest of this work. 



4 Results 

For Xi(rj) we have used the standard Kawano code (Kawano 1992) with N v = 3, r n = 889 sec and 
the correction AYp = +0.0006 (Kernan 1993). We begin by reviewing previous work on the gamma 
and log normal distributions and provide an extended analysis of the gaussian distribution. Then 
we consider a model based on inclusion of the scale of the perturbations where the power spectrum 
is given by a power law. 

4.1 Gamma Distribution 

The gamma distribution (EP, Yang et al. 1984) is given by 

f( v ) = v a - 1 e- a ^. (13) 

For this distribution the variance J 2 , is 



5r]\ 2 _ (rf) - 



, ' 2 =a-\ (14) 

77 J Tj 



The results of varying 5 2 and f) are shown in figure [l|. We have required the averaged abundances 
to fit the observations (||). Figure |l| shows the abundance contours of the light elements as given 
by the limits in (|5|) and thus delineates the resulting parameter space that reproduces the correct 
abundances. In Yang et al. (1984) the parameters 5 2 and f\ were constrained to fj ~ 3 and 5 2 < 3 
without using the 7 Li bound and a weaker upper limit on 4 He of Yp < 0.25. (For 5 2 < 1, the upper 
bound on f\ is relaxed to the homogeneous upper bound.) Here, as one can see from the figure, the 
more restrictive 4 He bound combined with the bound from D + 3 He yet with the weaker bound 
from 7 Li (||) allows us to constrain 2.8 < 7710 ^ 3.6. Including the WSSOK 7 Li bound this range is 
further constrained to 2.8 < 7710 < 3.3. These results are summarized in table |]. 

4.2 Log-Normal Distribution 

The log-normal distribution (Barrow & Morgan 1983) is given by 
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For this distribution 

fj = e^ W ' 2 . (16) 

The results of our search of parameter space are shown in figure Using the weaker 7 Li bound (Q) 
then fjiQ < 3.6 is allowed. This result is similar to that given by Barrow and Morgan (1983), though 
we are using more restrictive bounds on D + 3 He and 4 He. If we include the WSSOK 'Li bound 
we are restricted to fj\o < 3.2. These results are also included in table |j. 



4.3 Gaussian Distribution 

The gaussian distribution we consider is given by 

f(v) = exp - 



(v - 

2a 2 



(17) 



A previous study (Sale and Mathews 1986) considered a one parameter gaussian distribution with 
H = 0. Using this two parameter distribution we find 



r] = H + \ — a 



- M 2 /2<T 2 



vr 1 + erf(^/ v / 2t7 )' 



Note that this expression requires 



rj > \ — a. 

IT 



(18) 



(19) 



Thus a region of parameter space is already restricted by the mathematics. The results of the 
parameter space search are shown in figure ||[ Using the weaker 7 Li bound (||) then t/io < 3.7. If 
we include the WSSOK 7 Li bound then 7710 < 3.4. These results are summarized in table |]. 



4.4 Power Law 

We assume a featureless power law for the power spectrum (Peebles 1987a,b) 
Thus A 2 (/c) oc k n+3 and the average abundance (FT^) is 



X, 



dr}Vf(v) x i(v) 



1 



V 



n+3 n 

2 




drj 77/(77) 



n+3 

77 \ 2 



(20) 



(21) 



In the limit of 7? max - > 00 (& max — > 00) this reduces to the simpler case where the scale of the 
perturbations was ignored @. Hence for 7y max ^> f) we expect the results to be independent of the 
power law index, n. For each of the previously studied distribution functions we have performed 
the average as given above (|2l]). The results are shown in figs. The spectral index, n = —0.5 
was chosen for illustrative purposes and since it is the preferred value in the PIB model. 

As an example of the results with a power law distribution consider the gamma distribution 
(fig. ||). For 5 2 = 0.1 (a) using the weaker 7 Li (Q) there is a small band in the allowed parameter 
space between the 4 He and D + 3 He bounds. The WSSOK 7 Li value is only marginally consistent 
with the D + 3 He bound leaving a very narrow allowed region in parameter space. For 5 2 = 0.2 
(b) the limit from the weaker 7 Li value roughly overlays the 4 He limit and the WSSOK 7 Li value 
is inconsistent with the D + 3 He limit. For 5 2 = 0.5 (c) 7 Li is always produced in greater amounts 



6 



than allowed according to the WSSOK limit. Finally for 5 2 = 1.0 (d) 7 Li is always overproduced 
and the D + 3 He is inconsistent with the 4 He limit (there is no allowed parameter space). 

From our expression for rj c (k) (|TT|) , if ry max = 5 x 10~ 10 then A m j n 70 kpc. Thus for this 
77max the spectrum of perturbations must be cutoff at the (comoving) scale of 70 kpc. No structure 
on scales smaller than this could be created from baryon isocurvature baryon perturbations. To 
allow smaller scale structure to form we must lower A m i n and hence raise r/ max . Thus we use as 
our criteria for determining valid regions in parameter space that all the light element bounds are 
satisfied and r/ max > 5 x 10~ 10 . From fig. ||(a) we see that requiring structure on scales less than 
70 kpc eliminates the horizontal region where the light element abundances are consistent with 
observations. From figs, ^a-d), we see that this leads to a limit of rjio < 6.0 with the weaker 7 Li 
abundance and fjio < 4.0 for the WSSOK 7 Li bound. Similar results hold for the log normal (fig. |5j) 
and gaussian (fig. |6|) distributions. These results are also summarized in table p]. Notice, if we 
allow an arbitrary cutoff at 3.5 x 10 -10 < ?7 max ^5x 10~ 10 we can get t) to be as large as desired. 
However, this requires introducing new physics to explain the origin of the cutoff at 70-80 kpc. 

We noted above that the results should be largely independent of the power law index, n. To 
verify this we have examined the gamma distribution for three values of n. The values chosen are 
n = —0.5 (fig. |4j), n = —2.5 (fig. |7]), and n = 2.5 (fig. ||). Comparing these figures we find that 
decreasing n slightly shifts the bounds to higher f). This effect is most pronounced for 5 2 = 1.0 
(part (d) of the figures). However, since this effect is small and shift all bound in approximately 
the same manner, the limits quoted in table Q are valid for all values of n. 

5 Conclusions 

As shown in table [l] the extra parameters used in these models of inhomogeneous BBN allow for 
only a slight increase in rj over SBBN. This fact is easy to understand. First ignore the 7 Li limits. 
Since the abundance of 4 He is a monotonically increasing function of rj and the abundances of 
D and 3 He are monotonically decreasing functions of r/, when we include regions of high r\ we 
are overproducing (underproducing) 4 He (D and 3 He) in the universe. The slight increase in the 
allowed value for r\ comes from the fact that 4 He is a slowly varying function of rj. When we include 
the Li bound it provides the tight upper bound on rj. This is due to the fact that the observed 
abundance lies in the trough of the predicted BBN production (see WSSOK). Thus any regions of 
high r\ greatly over produce 7 Li and we cannot allow such regions to have a significant contribution 
in the universe. The slight increase in r\ allowed is due to the generous limits we have allowed for 
the Li abundance. 

The case of a power law distribution allows only a slightly extended range at the expense of 
adding two new parameters, n and r/ max - The reason for this again traces back to the above 
discussion of mixing in regions with too much or too little of the light elements. Furthermore, 
these results are essentially independent of n. This in turn places tight constraints on models that 
contain these perturbations. In models that allow the power spectrum to extend to small scales 
the density of baryons is restricted to the same region as homogeneous BBN. Models that impose 
a small scale cutoff can, if the cutoff falls in just the right region, lead to much higher values for the 
baryon density at the expense of adding new physics. Recently Gnedin, Ostriker, & Rees (1994) 
have found some non-gaussian baryon isocurvature models can circumvent these bounds. We will 
return to this topic in a future work. 
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Figures 



1. Parameter space plot for the gamma distribution. The acceptable parameter space is between 
the 4 He lines (solid), below and to the right of the D line (short-dash), above and to the right 
of the D + 3 He line (dash-dot), and below the WSSOK 7 Li line (long-dash) or below the 
Thorburn 7 Li line (long-dash, short-dash). The light gray shaded region satisfies the 4 He, D, 
D + 3 He, and Thorburn (weaker) 'Li bounds. The dark gray shaded region satisfies the 4 He, 
D, D + 3 He, and WSSOK 7 Li bounds. 

2. Parameter space plot for the log normal distribution. The acceptable parameter space is as 
defined in figure |]. 

3. Parameter space plot for the gaussian distribution. The acceptable parameter space is as 
defined in figure |l|. The region above the slanted solid line is not accessible for mathematical 
reasons (see text for details). 

4. Parameter space plot for the gamma distribution with scale for n = —1/2. The acceptable 
parameter space is as defined in figure |l[ In a) with 5 2 = 0.1, fj < 6 x 10~ 10 for r/ max > 
5 x 10 -10 . In b) with 5 2 = 0.2, the WSSOK 7 Li limit does not allow any region of concordance. 
In c) with 5 2 = 0.5, the Thorburn 7 Li limit falls below the D + 3 He limit. In d) with S 2 = 1.0 
7 Li is overproduced for all values of fj and r/ max and the region defined by the D and D + 3 He 
limits does not overlap the region defined by the Yp limits. 

5. Parameter space plot for the log normal distribution with scale for n = —1/2. The acceptable 
parameter space is as defined in figure |J. 

6. Parameter space plot for the gaussian distribution with scale for n = —1/2. The acceptable 
parameter space is as defined in figure [TJ. 

7. Parameter space plot for the gamma distribution with scale for n = —2.5. The acceptable 
parameter space is as defined in figure |]. 

8. Parameter space plot for the gamma distribution with scale for n = 2.5. The acceptable 
parameter space is as defined in figure 0. 
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Distribution 
Function 


Upper limits 


on 7710 


w/o WSSOK Y Li with WSSOK Y Li 


Gamma 


3.6 


3.3 


Log normal 


3.6 


3.2 


Gaussian 


3.7 


3.4 


Power law + gamma 


6.0 


4.0 


Power law + log normal 


6.0 


4.0 


Power law + gaussian 


5.0 


4.0 



Table 1: Upper limits on 7710 for various distribution functions (see text for details). 
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